Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN M?eX style file vl.4) 



On the Electric Field Screening by Electron-Positron Pairs 
in the Pulsar Magnetosphere II 

S. Shibata^, J. Miyazaki^, and F. Takahara^ 

^ Department of Physics, Yamagata University, 1-4-12 Kojirakawa, Yamagata 990, Japan 

^Department of Physics, Tokyo Metropolitan University, Minamiohsawa 1-1, Hachioji, Tokyo 192-03, Japan 

^Department of Earth and Space Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan 



(N 
O 

o 



0^ 



> 

O 
(N 
O 

o 



Accepted — ; in draft form 2001 September 14; revised on 27 May 2002 



ABSTRACT 

We study on the self-consistency of the pulsar polar cap model, i.e., the problem of 
whether the field-aligned electric field is screened by electron-positron pairs that are 
injected beyond the pair production front. We solve the one-dimensional Poisson equa- 
tion along a magnetic field line, both analytically and numerically, for a given current 
density incorporating effects of returning positrons, and we obtain the conditions for 
the electric-field screening. The formula which we obtained gives the screening dis- 
tance and the return flux for given primary current density, field geometry and pair 
creation rate at the pair production front. 

If the geometrical screening is not possible, for instance, on field lines with a 
super-Goldreich- Julian current, then the electric field at the pair production front is 
constrained to be fairly small in comparison with values expected typically by the 
conventional polar cap models. This is because (1) positive space charge by pair po- 
larization is limited to a small value, and (2) returning of positrons leave pair elec- 
trons behind. A previous belief that pair creation with a pair density higher than the 
Goldreich- Julian density immediately screens out the electric field is unjustified at 
least for for a super Goldreich-Julian current density. We suggest some possibilities to 
resolve this difficulty 

Key words: magnetic fields-plasmas-pulsars:general-gamma-rays:theory. 



1 INTRODUCTION 

One of the promising models for the pulsar action is the po- 
lar cap model where the field-aligned electric field acceler- 
ates charged particles up to TeV energies, and resultant cur- 
vature radiation generates an electromagnetic shower. The 
model can provide an environment for radio emission mecha- 
nism with a particle beam in electron-positron pair plasmas 
and may also explain the gamma-ray pulsation. The pairs 
will be a particle source of the pulsar wind. The polar cap 
potential drop is a part of the electromotive force produced 
by the rotating magnetic neutron star. The voltage is a few 
percent of the available voltage for young pulsars, while it 
becomes some important fraction for older pulsars. The lo- 
calized potential drop is maintained by a pair of anode and 
cathode regions, the formation mechanism of which is the 
long-standing issue of the polar cap accelerator. 

Because any deviation of the space charge density from 
the Goldreich-Julian (GJ) density poj causes the field- 
aligned electric field within the light cylinder, how the GJ 
charge density changes along a given magnetic field line 
plays an essential role in formation of the field-aligned elec- 



tric field. Note here that the GJ density is determined by 
the magnetic field geometry: pcj ~ —Q.B^ /2-nc, where Q. is 
the angular velocity of the star and is the field strength 
along the rotation axis. For instance, the outer gap, another 
promising accelerator model, is thought to appear around 
the null surface on which the GJ density changes its sign. 
Due to the field geometry, the outer gap originally has a 
pair of the anode and cathode regions without current and 
pair plasma. The space-charge-limited fiow for the polar cap 
model on field lines curving toward the rotation axis also 
has a similar GJ density distribution providing a pair of the 
anode and cathode regions (actually it is regarded as a ver- 
sion of outer gap with external current [Hirotani & Shibata 
2001]). 

Once pairs are created in the accelerator, since the pair 
number density far exceeds the GJ value |pGj/e|, where — e 
is the electronic charge, dynamics of the created pairs makes 
crucial effects in localizing the field-aligned electric field. 

For the polar cap models, the space-charge-limited flow 
is found to produce the field-aligned electric field in some 
cases where the current density and magnetic field geometry 
are assumed (Fawley, Arons & Scharlemann 1977, Scharle- 
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mann, Arons & Fawley 1978). Later on, effects of pair cre- 
ation and general relativity arc taken into account (Arons 
& Scharlemann 1979, Muslimov & Tsygan 1992). In these 
models, given a current density, the space charge density 
deviates from the GJ charge density, and as a result the 
field-aligned electric field develops. Since the work function 
of charged particles on the stellar surface is sufficiently small, 
either electrons or ions can freely escape from the surface, 
and therefore the electric field on the surface is assumed to 
be zero. Based on this scheme, primary electrons arc shown 
to be accelerated to energies high enough to induce an elec- 
tromagnetic shower. 

The current density on each magnetic linos of force will 
not be determined by the local electrodynamics in the po- 
lar cap, but will be determined in a global manner. Most of 
the rotation power of the pulsar is conveyed by the pulsar 
wind, which is accelerated near and beyond the light cylin- 
der. Although the acceleration mechanism of the wind is still 
mystery, it is certain that an energetically dominant process 
which determines the current distribution takes place near 
and beyond the light cylinder. When the wind is accelerated 
(or heated), the current runs across the field lines {E-j > 0) 
as a simple result of the Poynting theorem, provided that 
the ideal-MHD condition is more or less valid. Therefore, 
some of the current coming up from the star crosses the 
magnetic field lines corresponding to the wind acceleration 
and returns back to the star. The rest of the current reaches 
large distances (e.g., nebula shock), and the same amount 
with opposite direction comes back to the star. In any case, 
the current closure is seriously controlled by the dynamics 
near and beyond the light cylinder. The current determined 
in such a way runs through the polar cap. Thus, the cur- 
rent density distribution across the polar-cap magnetic flux 
should be linked to the outer magnetosphere: the current 
density distribution is most likely to be determined regard- 
less of the field geometry of the inner magnetosphere. Wo, 
therefore, treat the current density as an adjustable free pa- 
rameter in our series of paper (Shibata 1991, 1995, 1997, 
Shibata, Miyazaki & Takahara 1998): the current density is 
a free parameter for the polar cap model, and it is adjusted 
when the local model is linked with a model of the outer 
magnetosphere. 

For a space-charge-limited flow, it is shown by Shibata 
(1997) that both toward and away geometries axe possi- 
ble to generate a large field-aligned electric field causing an 
electron-positron pair cascade as long as the current density 
is reasonably large, comparable to the GJ current density 
Pgjc. In his paper, the cases of super-GJ current densities 
arc examined as well as the sub-G.l cases. Note that the 
GJ current density is simply the GJ current, /qj ~ /ifi'^/c 
divided by the polar cap area, where n is the magnetic mo- 
ment of the star. The GJ current is derived by an order- 
of-magnitude argument, and the net current circulating in 
the magnetosphere will not be far different from this value. 
However, if the current is not distributed more or less uni- 
formly over the polar cap, but is focused or fragmented in 
the magnetic flux tubes by analogy of planetary aurori, then 
the local current densities can be much higher than the GJ 
current density. Furthermore, the current density will change 
in a wide rage held lines by field line. 

These points motivate us to examine various combinar 
tions of current density and magnetic field geometry for the 



local accelerator. Although deviation of the space charge is 
found to cause the accelerating electric field, it is not clear 
how the electric field is screened to complete the localized 
potential drop. Only the known case in which a finite poten- 
tial drop is formed without pairs is the geometrical screen- 
ing, i.e., the case of toward geometry or general relativistic 
geometry (Muslimov & Tsygan 1992) with a certain current 
density which is sub-GJ. For other cases, there remains in 
general an unscreened electric field in a region where pairs 
arc created. Note that pair plasmas arc copious and pair 
polarization may become efficient. In this paper, we derive 
the condition for the electric-field screening immediately be- 
hind the pair creation front where unscreened electric field is 
assumed to exist, in the presence of returning pair. The con- 
dition can be used for any current density (super- or sub-GJ) 
and for both away and toward field line curvatures. We will 
apply the condition to the cases for which geometrical effect 
is not efficient, i.e., away curvature and supcr-GJ current 
density. Another important application is the case where 
pairs arc created before the geometrical screening completes. 

In the previous paper (Shibata, Miyazaki & Takahara 
1998, Paper I), we demonstrated that the pair polarization 
does work when enough pairs are injected, where we as- 
sumed that pairs are created at a single point; the required 
multiplication factor is ~ 10^ which is the same order of 
the value inferred from observations of the Grab pulsar and 
nebula. However, we suggested that the required multiplic- 
ity has to bo reached in a very small distance, and this will 
not be realized in an actual pulsar magnetosphere. 

We have made some simplifications in Paper I in order 
to make the physics as clear as possible. They axe, 1) pairs 
have the same initial Lorentz factor and outward momen- 
tum, 2) the pair creation occurs at a single point, and 3) 
no positrons return to the neutron star. In this paper, we 
relax the above assumptions 2) and 3) to study more realis- 
tic cases: electron-positron pairs arc produced continuously 
along field lines for a finite length, and some particles may 
stop before reaching the screening point and return to the 
neutron star surface. As in Paper I, we assume that no fric- 
tional force works between the various components. Under 
these assumptions, we calculate the electric field structure of 
the pair creation zone and determine the maximum electric 
field which can be screened under a presence of returning 
particles. Since the thickness of this zone is geometrically 
thin, we adopt one-dimensional approximation in this pa- 
per. 

In section 2, we formulate the problem, and in section 3 
we solve the one-dimensional Poisson equation analytically 
assuming that the pair creation rate is constant as a func- 
tion of the electrostatic potential. In section 4, we solve the 
problem numerically when the pair creation rate is spatially 
uniform to confirm that the analytic results of section 3 are 
sufficiently accurate. In section 5 we summarize our results 
and suggest a possible way out to resolve the screening prob- 
lem. 



2 MODEL AND BASIC EQUATIONS 

Let us consider a steady model of a field-aligned accelera- 
tor which has a finite potential drop with the electric field 
screened at the both ends. For definite sign of charge, let us 
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assume electrons are accelerated outwards, i.e., the electric 
field points toward the star. Gamma-rays emitted by the 
electrons convert into pairs beyond a certain surface called 
the pair production front (PPF); beyond PPF pairs are as- 
sumed to be created continuously in space. 

Pair polarization is an effect that screens out the elec- 
tric field. The pair positrons are decelerated soon after their 
birth, while the pair electrons are accelerated, and as a re- 
sult of continuity, a positive space charge appears to reduce 
the electric field. In this paper, we intend to evaluate the 
space charge produced by pair polarization in the presence 
of returning positrons. 

Notations are the same as Paper I: the strength of 
the magnetic field B is denoted by B, and the compo- 
nent along the rotation axis is B^; the non-corotational 
electric potential <E> gives the field-aligned electric field by 
i?!! = —B ■ V^/B; v and p are the velocity and density 
with the subscripts '1', and ' — ' indicating the primary 
particles, pair positrons and pair electrons, respectively; 
PG.T ~ — (r2_Bz)/27rc is the GJ charge density. We also use 
the following normalized quantities: 

e<I> 



potential 



effective 
space ctiarge 



Figure 1. An accelerating region with a finite potential drop has 
a pair of anode and cathode regions. For a parallel rotator (not 
anti-parallel), negative space charge appears in the side of the 
stellar surface, and positive space charge is required for 'screen- 
ing' the electric field at the other side of the potential drop. The 
positive screening charge can be provided by toward curvature 
(on field lines curving toward the rotation axis) and by a pair 
plasma. 
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where J = —enivi is the current density of the primary 
electrons; s is the coordinate along the magnetic field, which 
is normalized by the relativistic Debye length. 

The Poisson equation for the screening region may be 

_ ^ 



dE 
dx 



+ jo + P++ p- 



(1) 



where E = —d(j>/dx is the normalized electric field along 
a given magnetic field, /3i is the normalized velocity of the 
primary particles. In the right hand side of (|^), the first 
term represents the negative space charge produced by the 
primary electronic current stream {—j/Pi < 0), and jo rep- 
resents — Pgj, which is positive. The difference, —j/Pi + jo, 
is the space charge that produces the accelerating electric 
field. The values p± are the normalized charge densities of 
the pair electrons and pair positrons, which appear beyond 
PPF. Note that j and jo take positive values in our sign 
convention. 

In a steady acceleration region along given magnetic 
field lines, 'an effective space charge' which is represented 
by the right hand side of (|l|) becomes positive or negative 
at each side of the accelerating region as shown in Figure 1. 
The negative charge region is located near the stellar surface 
and electrons are accelerated upwards. In the space-charge- 
limited fiow (Fawley, Arons & Scharlemann, 1977) , the neg- 
ative charge is provided by non-relativistic electrons leaving 
the stellar surface; —j/(3i + jo < with /3i <^ 1. Additional 
negative charge is induced on the side wall of the fiow tube. 



As the electrons are accelerated, /3i 1, the effective space 
charge becomes positive above a certain height, —j+jo > 0, 
if j is smaller than jo- Thus, the anode region is formed in 
this case. 

Further development is to include the effect of magnetic 
field geometry: on field lines curving toward the rotation 
axis, jo = Bz/B (< 1) increases along the field lines while 
j oc J/B is constant, and therefore even if —j + jo < near 
the star, as the fiow goes out, — j -I- jo can be positive above 
some distance as far as j < 1 (sub-GJ). This toward curva- 
ture is efficient to screen the electric field which is sufficiently 
strong to cause pair cascade. 

One can choose a particular value of j so that the elec- 
tric field is localized with the screened electric field at the 
both ends without pairs. However, as has been mentioned, 
the current density is determined in a global way. It is there- 
fore unlikely that such a particular value of the current den- 
sity is realized in general. One can expect that PPF is formed 
at the place where the electric field remains unscreened. 

If j > 1, the anode formation by the geometrical effect 
does not work because — j+jo < 0. The same situation takes 
place for any value of j on field lines curving away from 
the rotation axis (Mestel 1981). In these cases, the steady 
accelerator is possible to exist only if an anode region is 
provided in some way for which we study the dynamics of 
pairs in the PPF region. 

Since screening by pairs occurs in a small scale such 
as 7mc^/e|-E|| I, j and jo are assumed to be constant in the 
screening region around PPF. 



3 CONDITION FOR SCREENING 

The screening region is schematically shown in Figure 2. 
The pair creation starts at PPF {x — 0), and continues 
beyond {x > 0). The non-dimensional electric field at a; = 
is denoted by ~E{, and the potential at PPF is used for 
the reference, i.e., S — nt x — 0. Let us assume that the 
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Figure 2. Schematic diagram of the screening region by pair po- 
larization. Electron-positron pair plasma is created at more or 
less uniform rate beyond the pair production front (PPF). The 
pair positrons are decelerated while the pair electrons are accel- 
erated. When the pair positrons become non-relativistic, there 
appears positive space charge, which should screen out the accel- 
erating electric field. Some of pair positrons will turn back to the 
star. Pairs created near the screening surface (SS) go through the 
region. 



field-aligned electric field is screened out at a certain point 
defining the screening surface (SS). 

Multiplying (0) by E and integrating it from PPF to 
SS, we obtain tfie screening condition 



2 ' 



(P+ + P- 



U -jo) 



(2) 



where we have assumed for the primary particles f3i — 1, 
and the subscripts 's' indicates values at SS. This condition 
is equivalent to that the column charge density between PPF 
and SS is equal to Ef (the Gauss' law). 

The net pair flux produced between PPF and a posi- 
tion X is given by the multiplication factor, m{x), times the 
primary electron flux, j. The multiplication factor m{x) is 
obviously a monotonic function of x. We assume that is 
also a monotonic function of x between PPF and SS. Thus 
either x, m or (p can be used to designate the coordinate 
position. 

The charge densities of pairs produced by a small flux 
element between m and m + dm at a 'position (j>' are given 
by 



dp± = ±jp = ± /^ dm, (3) 

where the Lorentz factors 7± are measured at the position 
(j) and are functions of the injection point denoted by ^(m), 
i.e.. 



for electrons and 

j+ ^ ~(p + (t){m) + 7 



(4) 



(5) 



for positrons; 7 is the Lorentz factor at creation, and as- 
sumed to be a universal constant, say 7 ~ 10^. 
Positrons return back to the star when 



(m) + 7 - 1 < 



(6) 



where (jirctim) indicates the returning 'position', otherwise 
positrons have enough energy to go through SS. Thus there 
is a separatrics Mi below which the positron flux (0 < m < 
Ml) returns while beyond which the positron flux (A/i < 
m < M2) passes through SS, where M2 is the multiplication 
factor at SS. 

The potential at which the flux Mi is created is given 
by the marginal condition of (P) 

,^1 =,^(Mi) = 0,-7 + 1. (7) 

If there are returning positrons, the space charge be- 
fore PPF becomes — (j//3i) + jo + jMi, so that positive 
contribution is added by the returning flux. For accelera- 
tion of the primary electrons, the column charge density 
between the star surface and PPF must be negative. If too 
many positrons are reflected back, the negative column den- 
sity below PPF will be broken down. If Mi is given and 
is small enough to assure the negativeness of the column 
charge below PPF, the growth of the electric fleld before 
PPF is calculated for an assumed j as in the conventional 
space-charge-limited flow. 

Once the primary particles become relativistic, the 
space charge before PPF is represented by —j + jo + jMi. 
If this space charge is negative, then we will have further 
acceleration, while if it is positive, then screening. We may 
have —j + jo + jM\ > 0, for instance, on field lines curving 
toward the rotation axis. The positiveness does not guaran- 
tee vanishing of E\\ at PPF because PPF may be formed at 
a place where E\\ still remains unscreened. For this case, the 
present analysis gives the screening distance after PPF. On 
the other hand, the PPF can be formed in a negative space 
charge region where 



- j + 30 + jA/i < 0. 



(8) 



This takes place on field lines with away curvature or with 
super-GJ current densities (for any type of field line curva- 
tures). For this case, screening relies on pair polarization. 

It is notable that the observed X-ray flux from the 
heated polar cap suggests that Mi is much less than unity. 

In addition to the screening condition (2), we have a 
supplementary condition that the second derivative of cj> 
must not be positive on SS: because d<j)/dx is positive at 
X — and zero at SS, if d^(j)/dx^ were positive at SS, dcp/dx 
would be negative just below SS, indicating d(l>/dx vanishes 
inside SS, which contradicts the definition of SS. Integration 
of (^ up to M2 yields the pair charge density at SS, so that 
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we have positiveness of the effective space charge density at 
SS as 



M2 



dm 



M2 



dm 



Ml 



(9) 



With the help of the contribution of pairs to 'col- 
umn charge' in the </)-coordinate is given by 



(P+ + P-)d<t) 







fMx l-<t>rct{m) 

jdm I 

10 J<l>(m) 

+ / jdm 
Jo Jo 

f-A/i 

jdm 

10 J t/>{m) 

/•A/2 

+ / jdm 

J Ml J,l>{m.) 



1 1 



d(j). 



(10) 



In the above equation, the first term comes from the up- 
going positrons which eventually stops at (j>rct, the second 
term is due to the returning positrons. The third term is for 
the electrons which are paired with the reflected positrons. 
The fourth term is for the pairs which pass through the 
region. 

The first term is integrated exactly to yield 



1st term 

0+ ^ 



7+ 



jdm / 

J4>(m) 



zdfj) 



jMivZ-y'-l 
jMi7-jA/iAi, 



(11) 



where the last expression is obtained if 7 ^ 1, and Ai 
(27)"^ <C 1. The second terms becomes 



2nd term 

^ 0+ 



Ml 



jdm^J + ct>{m)Y 



rMi 
J I dm 



(7-h0(m))- 



2(7 + 0(m)) 



Ml 



/•Ml 

-j / [4>s - (f>{m)]dm + jMi A3 (14) 
Jo 



with 



1 \"'^ 1 



(15) 



where 7_s(m) = j + <j)s — <j}(m) is the Lorentz factor of pair 
electrons at SS. 

The fourth term becomes 

= j dm 

J Ml 



4th term 



2^72 _ 1 „ ^j_2^ _ 1 



(16) 



where, 7+s(m) = 7 — 0s + (l>{m) is the Lorentz factor of the 
positrons at SS. If one assumed 7, 7-a and 7+a S> 1, then the 
integrand would be 7— 7-s+7— 7+s which would vanish after 
inserting the definition of 7_s and 7+5. This means simply 
that pairs passing through with relativistic speeds would not 
cause any space charge. However, the positron flux near Mi 
causes positive net charge because 7+s ~ 1 near </!> 7- (j>s. so 
that the 4th term cannot be ignored: we rearrange (tq) as 

M2 



4th term 



l\/ 7^ — 1 — 7]dm 



Ml 

M2 



7+s]dm 



Ml 

M2 



1 - 7_a]dm, (17) 



Ml 



which can be integrated exactly (see Appendix) if the pair 
creation rate is constant with respect to the '^-coordinate'. 
In this approximation, we have 



m = Ml + [4>{m) — ^i] 



with 



M2 - Ml M2 - Ml 



7- 



1 



= const. 



(18) 



(19) 



jMi^ + j j (j)(m)dm- j Ml A2 (12) 
'0 



where is the pair multiplication factor created in a unit 
potential drop m(? je, in other words, is the pair mul- 
tiplication factor produced in the 'braking distance' As = 
mc? le\E\\\. After integration we have 



with 

A2 = 



2(7 + <^(m)) 



« 1, 



(13) 



4th term , . ^ 



(20) 



where {/(?ti))™J indicates an average with respect to the 
flux m, defined as 



7712 — mi 



dmf{m) 



The third term of (|l^) is the column charge by electrons 
associated with the returning positrons: 



3rd term 
e~ — > 



l-Mi 

-J I dm 



V(7 + 0s-0(m))2 



V7^-l 



where ^ is a constant of order of unity; ( — 1.744 for 7 = 100 
(see Appendix for derivation) . 

Near SS, d(f>/dx — > 0, = dm/dcj) is not constant if 
the pair creation is spatially uniform. In the next section, 
we perform a numerical calculation for a spatially uniform 
pair creation and find that the above evaluation is quite 
accurate. The analytical evaluation by using a constant rate 
acf, = dm/d(j> will be appropriate for more or less uniform 
pair creation. 

Inserting the above evaluation of all the four terms into 
(^) and rearranging relevant terms, we have 

^Ef = 2jMi - 2jMi {AcP) + ja^C - (j - jo - jMi)^, (21) 
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where 



Ml 



{A<j)) = — / [(j)s - (l>vct{m)]dm 



^ Jh J 0(m)]dm 



(22) 



is the potential drop (flux weighted) in the region where 
returning lasts, in other words, the potential difference be- 
tween where returning start and where screening completes. 

The pair effect is made up of three components corre- 
sponding to the first three terms in (|2l|). The first effect is 
a positive charge produced by non-relativistic part of the 
returning positrons. The second effect is negative because 
it is due to electrons left behind the returning positrons. 
This effect by returning is fatal: the more positrons return, 
the more electrons are left behind to create negative space 
charge in the downstream (outer) region, whereas the re- 
gion needs positive charge for screening. The last effect is 
caused by polarization of the pairs which pass through SS 
(3rd term). 

If the region of returning is large, in other words, if 
(Atj)) ^ 1, the pair screening becomes obviously difficult 
to complete because the negative space charge produced by 
pair electrons left behind in upstream region becomes large 
to prevent screening. If uniform pair creation is assumed, 
(Af6) = (l/2)(^i, and Mi = a^cj)-i_. In this case, the condition 
(Ell) becomes 



\Ej = ja^ [01 (2 - 0i) + CI - (j - JO - jMi )(^, 



(23) 



The first term of the right hand side of (^3|) corresponds to 
the pair polarization. The second term is the 'background 
charge', the sign of which depends on the global parameters 
j and jo, and where PPF is formed. Let us restrict ourself 
to the case for which the global effects is not helpful but 
pairs cause screening, i.e., the last term makes no positive 
contribution. In order for the right hand side of ( pSj ) to be 
positive, the necessary condition is that the pair polarization 
term is positive, and then we have 



2(A<?i) = <^i < 1 + ^Jl + C ~ 2.6, 



(24) 



which means that returning of positrons lasts only a few 
braking distance. 

Here we invoke the condition (^ that is for charge den- 
sity at SS to be positive: 

. /"^^^ dm . f^^^ dm , . 

Ps = -J ^ + J / -3 (j - Jo) 

Jo P- Jaii P+ 

= ja^ (2V7' - 1 - v/(27 + 'i^'i - 1)' - 1) - (j - io) 

« (1 - 0i)a0j - (j - Jo) 

= ja^(l-20i)-(j-jo-Mij) >0. (25) 

As was mentioned, returning positrons leave negative space 
charge at SS, the positiveness of pa sets an strong constraint: 
from (^5|), again noting no positive contribution from the 
last term, we have (1 — 2(j>\) > 0, i.e.. 



< 01 < -. 



(26) 



[j - jo- jMi)4>, 



(27) 



where C,' is between C, = 1.744 (when 0i = 0) and C + 3/4 = 
2.494 (when 0i — 1/2). Eventually the necessary condition 
for pair polarization to screen out E{ becomes 

> (28) 

Recall that is the pair multiplication factor produced 
within the braking distance As = mc^ /e\E^^\. In dimensional 
form, the required pair multiplication factor in As is given 
by 

> J^^. (29) 



STrmc^noC'j 



For polar cap models, typical values of the right hand side 
of (pS) are as large as 10^. For instance, 
acceleration of 7 ~ 10^ in a distance of L • 



given a primary 
10^ cm, we may 



have l^iil = 1.7 x 10^ in esu or K/Stt = 1.2 x 10^(76/1-4) 



erg cm 

AMsnrcr 



at PPF, and 



>2.0xl0^(f )^(9Li)^ 



(30) 



where 75 = 7/IO® , and L4, Po.i and Bu are in units 
of 10'* cm, 0.1 sec and lO^^G, respectively. This multi- 
plication factor should be achieved within a distance of 
As = m(? / e\E\f^ \ = 7~^L ~ lO"'^ cm for screening. It seems 
quite difficult that the pair multiplication factor of order 
of 10^ in this sort of small distance. If the pair production 
takes place over a dimension of say 10* — 10^ cm, then the re- 
quired multiplication factor in the whole region is 10^ — 10*" 
and is much higher than the value predicted by conventional 
magnetic pair production. If, on the other hand, we assume 
L ~ 10® cm, then AMscrocn ~ 0.2, As = 1 cm and the total 
multiplication factor becomes 2000 to 20000, which may be 
allowed. 

In the case where j < jo, the screening completes with 
the help of the 'background' positive space charge, i.e., — (j — 
jo — jMi)(f>s in (p^). However, the distance to the screening 
region is enlarged by the negative contribution due to the 
pair electrons left behind. For a given PPF, i.e., given j, jo 
and the pair creation rate a^, we can calculate the screening 
distance in potential scale and return fiux Mi in use of 
(H), (|) and Afi =a^0i: 

^ _ il/2)Ef+ja^if- 



■C) 



(31) 



The condition (24) is then rewritten as 



jo - j[l - a0(7-h 1)] 

It can be seen that the screening is possible if j < jo /[I — 
a^{"f + 1)] = jcr- For the conventional polar cap model, 
j^^ jo because a^^ is much less than unity. If j — > jcr. Mi 
becomes large so that we have to treat the whole acceleration 
region including the cathode region, which will be discussed 
in a subsequent paper. 



4 NUMERICAL RESULTS 

In the previous section, we have discussed the screening con- 
dition under the assumption that the pair creation rate is 
constant with respect to the electrostatic potential. In this 
section, we present numerical solutions when the pair cre- 
ation rate is spatially uniform, dm/dx = constant. 
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Figure 3. Electrostatic potential as a function of the distance 
from the pair production front. 
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Figure 4. Electric field strength as a function of the distance 
from the pair production front. 



First we check the code by calculating the same case as 
in the previous section. We present the case of 7 = 100, j — 
2jo, Q0 = 1 and Mi — 0.5, and hence (f>i — 0.5. The results 
are shown in Figures 3 and 4, where <j> and \E\ as functions 
of X are shown. We find that the screening of the electric 
field occurs in the region where non-relativistic positrons 
appear, and the distance scale of this region is roughly the 
braking scale of mc^ /e\E\\\. The whole dimension between 
PPF and SS is 7 in potential scale. Consequently 99% of 
the pairs produced in this region go through and only 1% of 
positrons return back to the star. In this case, the numerical 
calculation gives E{ = 3.159, which perfectly agrees with the 
numerical estimate based on (|2^), 

Next, we examine the case where the pair creation rate 
is spatially constant, i.e., m{x) — UxX with a constant ax- 
For 7 = 100, i — 2 jo, and Mi = 0.5, we solve (Q) numeri- 
cally. When — 0.5, we obtain the electric field strength 
at PPF, Ef — 3.160, which is almost the same value as the 
case of the constant a^. If one plots as in the same diagram 



120 




Figure 5. Multiplication factor as a function of 4> when dm/dx 
is constant. 

as Figures 3 and 4, it is hard to find the difference on the 
plots. 

The small difference found in Ei between the two cases 
arises from the behavior of the pair creation rate near = 
08. When we assume m = a^<j), the spatial pair creation 
rate becomes slightly smaller as the electric field becomes 
weaker near SS. With this numerical solution, we find that 
( p8[ ) is useful even for the case of the spatially uniform pair 
creation. 



5 CONCLUDING REMARKS 

In this paper, we considered the electric field screening by 
polarization of electron-positron pairs which are created be- 
yond the pair production front as is generally supposed for 
the polar cap models. We allow the pair positrons to return 
back to the star as far as the condition for the accelera- 
tion is satisfied. The one-dimensional Poisson equation is 
solved together with particle motion for a steady state ac- 
celeration region. We obtain the formula (|2^) which give 
the electric field strength that can be screened for a given 
current density j, a given field geometry jo and a given pair 
creation rate a^. It is found that (1) pair polarization has lit- 
tle contribution for positive space charge, and (2) returning 
of positrons makes the screening difficult seriously because 
the pair electrons left behind the returning positrons pro- 
duce negative space charge in the screening region where the 
positive space charge is required. As a result, the thickness 
of the screening is restricted to be as small as the braking 
distance As ~ mc^/e|_Ej|j for which positrons become non- 
relativistic. We confirmed the previous result of Paper I that 
the electric field in the case where geometrical screening is 
not possible, i.e., —j + jp + jMi < at the pair production 
front, we have, from (ES), 

i < (^) ^'^-^^-^ (^2) 

where AMscrccn is the pair multiplication factor within As. 
If the primary current density is of order of the Goldreich- 
Julian (GJ) value, the required pair multiplication factor 
per one primary electron is enormously large and cannot be 
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realized in the conventional pair creation models. A previ- 
ous belief that pair creation with a pair density higher than 
the GJ density immediately screens out the electric field is 
unjustified. 

Some mechanism to salvage this difliculty should be 
found. As already mentioned, if the last term of (p3|), 
— (j — Jo — jMi)(j}s, becomes positive screening is possible 
as far as j < jcr ~ jo- Since jo — B^/ B < 1, j must be less 
than unity. It must be noted that the positive space charge 
exists before the pair production front in this case, so that 
screening takes place before the pair production front and 
is reduced. Even a reversal of is possible for a small j 
(Shibata 1995). In any case, if geometrical screening works, 
a finite potential drop along field lines can complete. But 
the current density is restricted in a certain range. 

The condition (|2^) implies that if the current density is 
much higher than the GJ value, the required multiplication 
factor can be reduced. Although the total current circulat- 
ing in the magnetosphere should be of order of so-called GJ 
current, 7gj = fiQ^ /c, where fi is the magnetic moment of 
the star, current densities can be in principle focused: the GJ 
current density is just a value if 7gj is more or less uniformly 
distributed over the whole polar cap area, and therefore if 
the current distributed, say, in an annular region with fine 
structures such as seen in aurora, the focused current density 
will be much higher than the GJ current density. Thus, if j is 
much larger than jo, and the required multiplication factor 
will be reduced. (If the current is super-GJ, toward curva- 
ture will not helpful for screening because — (j — jo ^jMi) is 
negative.) In any case, physics in focused field-aligned cur- 
rent, e.g., efi^ects of inhomogeneity across the magnetic field 
lines, two stream instability in a high current density with 
unscreened electric field, has yet been properly studied. 

Since 7 is an order of 100, (f>s turns out to be an order 
of 100 using (0) and (H). From ^ we may safely approx- 
imate as = A/2/100. From (|2^, E{ < \/M2/5, taking 
— 2 and j = 1. As an example, let us assume that total 
pair production occurs with a multiplication factor 10* over 
the distance Xp = W^. Thus, M2 = 10''(a;s/a:p) = Xb- Noting 
that E{ « (jjs/^s = 100/a;s, we obtain Xs > 60 and Ef < 1.5. 
Since the unit length c/ujp is typically 1 cm, this electric 
field corresponds to a voltage mc? /e per cm. Assuming that 
primary electrons are accelerated up to 7 > 10^ with this 
weak field Ef ^ 1, the required length between the neutron 
star surface and the pair production front is about 10^ cm in 
hight. Thus pair screening is possible if the pair production 
front is located high above the surface and if the field-aligned 
electric field is not so strong. However, the previous models 
predict a stronger electric field, so we would need an addi- 
tional ingredient to keep a small field as long as a stellar 
radius. In any case, a weak electric field kept in a long dis- 
tance is one possible solution to have a self-consistent polar 
cap model. 

Another possible way out is to include frictional forces 
between various components of charged currents. Since fric- 
tion on positrons pulls them outwards along with electrons, 
returning fraction of positrons will be reduced, which will 
make screening easier. Although there are some ideas for 
physical processes of friction (two stream instability, pro- 
duction of positroniums and others) , it is not clear whether 
the frictional force can be strong enough to lift positrons and 



screen the electric field. Study on the frictional interaction 
under unscreened electric field is strongly demanded. 
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APPENDIX A: SPACE CHARGE BY PAIRS GO 
THROUGH THE SCREENING REGION 

The space charge by pairs that go through the screening 
region is given by (p^). The second and third terms can be 
integrated analytically if the pair flux is proportional to the 
potential for which 



dm — a^d(j){m) 



(1) 



where the pair creation rate in unit potential drop is 
assumed to be constant. 

The second term of ( |l7| ) becomes 

j J dm (^7+s - 



(2) 



J<t>i ^ ' 
= - 7V^7^ - 1 + ln(7 + V^7^ - 1) - 1 

In the same way, the third term of ( [1^ ) becomes 

j I (7-s - V7-S - 1) a^d(j>{m) 

-7-s/^E~T 



f A/2 

J I dm I 

' Ml 

r4'^ 



+ ln(7-, + ^72^-l)-l 



27-1 



(3) 



Thus we arrive at (|^) with ( = 1.744 if 7 = 100 is 
substituted. 
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